202212.00150v1 


chinaXiv 


第 39 卷 第 5 期 
2022 年 9 月 


干枯 骆 怠 刺 对 风沙 流 场 


XE, 


FRE, 


干旱 区 研究 
ARID ZONE RESEARCH 


FRAS 


ChinaXiv 合 作 期 刊 


Vol.39 No.5 
Sep. 2022 


影响 的 数值 模拟 研究 


KER’, 


aR i 


(1. 北京 林业 大 学 水 土 保持 学 院 , 北 京 ” 100083; 2. 新 疆 维吾尔 自治 区 水 土 保持 生态 环境 监测 总 站 ， 
新 疆 乌鲁木齐 830002; 3. 喀什 大 学 生命 与 地 理科 学 学 院 , 新 疆 喀什 844000) 


摘 要 : 植物 固沙 是 荒漠 地 区 防治 风沙 灾害 的 重要 措施 之 一 ,骆驼 刺 (41pagi camelorum) 作 为 典型 荒漠 植物 具有 重 
要 的 应 用 价值 。 利 用 fluent 软 件 对 30 cm 高 度 干 村 骆驼 刺 附 近 的 流 场 进行 数值 模拟 ,分 析 风 速 特征 和 积 沙特 征 , 并 
加 以 野外 试验 验证 。 结 果 表 明 :(1) 植株 附近 流 场 大 致 可 分 为 遇 阻 减速 区 .抬升 加 速 区 . 色 流 减速 区 以 及 恢复 区 。 


植株 后 会 形成 微弱 的 涡流 ,涡流 的 回流 区 高 度 与 距 植株 远近 有 关 ,但 整体 低 于 0.14 mo (2) 当 风 速 为 6m's 时 ,植株 
主要 影响 0.6 m 高 度 以 下 的 水 平 风 速 。 植 株 后 一 定 距离 内 的 水 平 风 速 随 高 度 增加 不 在 呈现 严格 的 对 数 分 布 ,而 是 
存在 2 个 极 小 值 ,并 在 0.3~0.6 m 高 度 范 围 内 以 较 大 的 加 速度 快速 增 大 ,(3) 植株 的 防风 效率 整体 呈现 出 随 风速 增 大 
而 减 小 的 规律 , 且 该 现象 随 高 度 增加 愈加 明显 。 当 风速 由 6 ms "' 增 大 至 10 ms "时, 株 后 5.3 m 范 围 内 ,0.3 m 高 度 
处 的 防风 效率 由 40% 减 小 至 16.56%。(4) 植株 附近 的 积 沙 因 风 速 的 不 同 具 有 差异 , 当 风 速 较 小 时 , 积 沙 主 要 集中 在 


前 方 植株 附近 以 及 植株 之 间 , 随 着 风速 的 增 大 , 积 沙 后 移 。 
>; 流 场 分 布 ; 数值 模拟 ; 骆驼 刺 ; 涡流 
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受 干旱 半 干 旱 气 候 的 影响 ,风沙 流 成 为 我 国 西 
部 荡 漠 地 区 近 地 表 物质 运动 的 主要 表现 形式 " ,其 
进一步 发 展 极 易 导致 干旱 半 干 旱地 区 生态 恶化 , 引 
起 风沙 灾害 ”, 在 这 一 灾害 中 ,地 面 沙 粒 构成 了 物质 


研究 ,例如 毛 泽 叙 ”利用 多 孔 介质 模型 ,对 森林 影响 
下 的 风 场 分 布 特征 进行 模拟 , 唐 朝 胜 等 "探讨 了 术 
胶 防护 林 的 防风 效应 ,这 一 类 研究 主要 是 针对 风声 
变化 进行 分 析 , 对 于 沙 粒 的 移动 与 沉降 关注 较 少 。 


基础 , 风 提 供 塑 造 地 表 的 动力 条 件 , 也 是 风沙 灾害 
的 直接 动力 来 源 “。 植 被 以 其 综合 优势 在 覆盖 地 
表 、 分 解 风 力 、 阻 挡 输 沙 等 多 方面 对 地 表土 壤 风 人 刨 
和 风沙 流 结 构 产 生 影 响 ” ,植物 固沙 技术 作为 一 种 
可 以 长 效 改善 豪 漠 化 的 措施 ,在 防风 固沙 技术 中 占 
据 主 导 地 位 。 

目前 ,对 植物 防风 固沙 效益 的 研究 主要 采用 3 
种 方法 ,一 是 传统 的 野外 观测 ,二 是 模拟 实验 ,包括 
风 洞 模拟 实验 和 数值 模拟 实验 ,三 是 基于 遥感 -地 
理 信息 系统 的 研究 。 数 值 模拟 仿真 技术 可 以 大 大 
减少 实验 工作 量 , 同 时 获取 野外 试验 难以 获取 的 数 
据 , 随 着 近 些 年 的 不 断 发 展 ,被 广泛 应 用 于 流体 力 
学 、 风 沙 物理 学 以 及 风沙 环境 等 多 种 领域 。 利 用 数 
值 模 拟 研 究 植 物 防风 固沙 主要 通过 2 种 模型 ,一 种 
是 多 孔 介 质 模 型 ,主要 应 用 于 和 森林、 防护 林 的 流 场 
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另 一 种 是 采用 实体 几何 直接 开 孔 建立 模型 ,多 应 用 
于 对 植株 "以 及 植物 沙 障 ( 草 帘 沙 障 . 草 方 格 沙 
障 .芦苇 沙 障 等 ) 防风 阻 沙 性 能 的 研究 ,例如 来 
风 兵 ”模拟 研究 了 单 株 胡杨 (Populus euphratica ) JE 
围 流 场 结构 特征 , 陈 柏 羽 等 根据 芦苇 (Phragmites 
australis ) 形 态 特征 ,建立 简化 模型 ,分 析 上 芦苇 沙 障 周 
围 风 速 变化 及 积 沙 分 布 特征 ,并 得 到 了 合理 布设 
间距 。 

W YERI (Alhagi camelorum) 是 一 种 典型 的 还 漠 
植物 ,其 根系 发 达 , 生 命 力 项 强 , 是 我 国 西部 充 漠 区 
最 常见 的 防风 固沙 植物 之 一 。 骆 驼 刺 以 单 株 或 者 
灌 从 的 形式 存在 ,不 同 区 域 骆 驼 刺 的 大 小 具有 差 
异 , 较 大 的 骆 弦 刺 灌 从 会 使 大 量 沙 粒 沉降 堆积 ,经 
过 长 时 间 的 风沙 活动 形成 灌 从 沙 堆 。 每 年 5 一 7 月 
是 骆驼 刺 的 生长 旺季 ,进入 生长 季 的 骆驼 刺 ,新 生 
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梳 芽 在 前 一 年 枯死 的 旧 术 根 部 重新 发 出 , 枢 叶 密度 
增加 ,高度 和 冠 幅 增 大 ,区 域 的 植被 覆盖 度 相应 增 
大 ,8 月 以 后 植株 的 木质 化 过 程 使 得 地 上 部 分 在 冬 
季 干 枯 后 , 除 末 端 细小 碎 枝 脱落 外 ,大 部 分 枝条 直 
立 不 倒 , 直立 的 植株 可 以 消耗 风 动 量 , 减 小 剪 切 力 ， 
虽然 春季 活体 骆驼 刺 植株 非常 有 限 ,但 村 立 的 植株 
仍然 可 以 产生 良好 的 防风 阻 沙 效果 。 根 据 野 外 
调查 ,本 文 研究 区 春季 干枯 骆驼 刺 平 均 高 度 为 30 
cm 左右 ,分 布 方式 为 零散 分 布 ,研究 表明 ,风沙 运动 
中 沙 粒 的 运输 主要 集中 在 近 地 表 30 cm 内 “”。 通 过 
数值 模拟 ,对 典型 充 漠 地带 春季 干枯 骆驼 刺 对 近 地 
表 处 风沙 活动 的 影响 进行 研究 ,分 析 植 株 附近 流 场 
的 变化 ,并 利用 野外 试验 对 模拟 结果 进行 验证 ,为 
利用 骆驼 刺 乃 至 植物 措施 防风 治 沙 提供 依据 。 


1 数据 与 方法 


1.1 研究 区 概况 

野外 试验 在 中 国 科 学 院 策 勒 沙漠 研究 站 实验 
观测 场 进行 。 策 勒 站 位 于 策 勒 县 ,在 塔里木 盆地 南 
缘 1400 km 风沙 线 上 。 策 勒 县 位 于 新 疆 西 南部 , 昆 
CAGE ,塔克拉玛干 沙漠 南 缘 ,地 理 位 置 为 33"18'~ 
39°30'N 、80°03'~82°10'E。 气 候 为 极端 干旱 性 大 陆 
匠 漠 气候 ,干燥 ,昼夜 温差 大 ,日 照 时 间 长 ;水 资源 
短缺 ,年 平均 降水 量 35.1 mm ,年 潜在 蒸发 量 25$95.3 
mm, 水 资源 补给 以 昆仑 山区 融雪 为 主 ;年 均 气 温 
11.9 %C ,极端 高 温 41.9 C ,极端 低温 -23.9 % ,无 霜 期 
209 d。 生 态 系统 脆弱 ,沙漠 、 苞 壁面 积 达 95% , AK 
植被 以 多 年 生 荡 漠 植 物 为 主 , 呈 斑 块 状 分 布 ,有 骆驼 
刺 、 花 花 柴 (Karelinia caspia ) FEA Tamarix chinensis) 
等 混和 后 群落 ,总 盖 度 小 于 15%""。 由 于 策 勒 处 于 塔 里 
木 盆 地 两 大 主导 风向 (NW NE) 的 下 风 区 域 ,全 年 盛 
行 西北 风 , 风 沙 活动 频繁 , 沙 侍 日 数 多 。 年 平均 风 
速 1.9 m's ,年 均 8 级 以 上 大 风 为 3~9 次 ,年 起 沙 风 
次 约 90 d ,土壤 质地 轻 且 易 形 成 风沙 流 ,沙尘暴 天 和 气 
高 发 ,最 多 可 达 59 d, 扬 沙 、 序 尘 天 数 达 240 do HF 
春季 气温 回升 快 , 冷 空气 活动 频繁 ,该 时 期 风沙 活 
动 尤 为 强烈 22 。 
1.2 数值 模拟 方法 
1.2.1 几何 建 模 和 网 格 划分 “基于 前 人 对 植被 以 及 
植物 沙 障 的 建 模 研究 ”” ,本 文 尝 试 建立 春季 干枯 
状态 骆驼 刺 的 简化 模型 ,由 于 多 孔 介 质 模型 参数 不 
易 设 定 , 因 此 采用 实体 开 孔 模式 ,直接 建立 孔 际 几 
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何 模型 。 风 沙 流 是 沙 粒 在 多 种 作用 力 下 的 三 维 运 
动 现象 ,根据 沙 粒 在 运动 过 程 中 主要 受 重力 和 拖 搜 
力 ,作用 力 基 本 处 于 同一 平面 上 ,因此 建立 二 维 简 
化 模型 。 综 合 考虑 野外 试验 和 本 文 的 研究 内 容 , 进 
行 合理 试 算 ,计算 流域 长 度 设置 为 20m, 高 度 为 3 
m, TE iL AE BOR s SER AY TB 30 cm, MAREE 
上 边界 2.7 mm, 按 照 形 态 进行 简化 ,经 多 次 建 模 模拟 ， 
对 比 计算 量 及 计算 结果 得 到 植株 模型 ,在 距离 入 口 
5m 和 6 m 处 各 放置 1 株 , 流 场 模型 简 图 如 图 1 所 示 。 


图 1 流 场 模型 简 图 
Fig. 1 Model diagram 


模型 计算 域 网 格 划 分 类 型 采用 四 边 形 网 格 
(Quad) ,划分 方法 采用 结构 化 网 格 划 分 。 由 于 风沙 
流 受 边界 层 的 影响 较 大 , 沙 粒 主要 集中 在 近 地 表 
处 ,因此 对 近 地 面 以 及 植株 附近 的 网 格 进 行 局 部 加 
密 。 近 地 面 划分 20 层 边界 层 网 格 ,植株 两 侧 划分 16 
层 边 界 层 网 格 , 第 一 层 网 格 尺寸 设置 为 0.015 mm, 
左右 以 1.05 的 比例 递增 ,计算 域 共 划分 网 格 数量 
为 19500 个 ,最 小 正 交 质量 (Minimum Orthogonal 
Quality) 远 远大 于 0.99, 最 大 正 交 焉 斜率 (Maximum 
Ortho Skew) 远 远 小 于 0.001, 网 格 质 量 良 好 ,满足 计 
算 要 求 (图 2)。 

12.2 边界 条 件 和 计算 参数 ”根据 流体 马赫 数 小 于 
0.3 时 为 不 可 压缩 流 的 空气 动力 学 原理 ,风沙 流 为 不 
可 压缩 流 。 计 算 域 左 侧 风速 进口 边界 条 件 为 速度 
A H (Velocity-Inlet) ,风速 计算 见 公式 (1); 由 于 出 口 


2 计算 区 域 网 格 划分 示意 图 


Fig.2 Schematic diagram of grid division of computing area 
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流速 和 压力 均 未 知 , 且 流动 完全 发 展 , 因 此 右 侧 边 
界 出 口 条 件 为 完全 发 展 自由 出 口 (Out-Flow) ; 对称 
面 上 不 存在 扩散 通 量 ,计算 域 上 壁面 边界 采用 对 称 
边界 (Symmetry) ,可 以 在 保证 精度 的 条 件 下 ,加 快 计 
算 速度 ;植株 模型 和 地 面 即 下 壁面 边界 采用 固体 无 
滑 移 边界 (Wall) ,粗糙 度 取 0.002 m, 

介质 沙 粒 为 稀 相 ,风沙 流 携 沙 粒 径 一 般 为 
0.075~0.25 mm, 将 沙 粒 近 似 为 单一 的 颗粒 球 模型 ， 
沙 粒 粒 径 ds=0.1 mm, 初 始 沙 粒 体积 分 数 为 0.02%， 
类 型 为 流体 (Fluid)”。 沙 粒 密 度 p,=2650 kg'm”, Rh 
度 凡 .=0.047 Pas。 空 气 密度 p=1.225 kg*m ,空气 运动 
Ft FE u=1.7894x10° Pas, 压力 为 标准 大 气压 ,重力 加 
速度 g=9.8 ms"。 入 口 边 界 速度 为 典型 风速 廊 线 流 : 

v(y)= nn (1) 

UE : vo Ay BE BRE 5 yo A HA Bit ERE 5k AG RIT A 
数 , 取 值 为 0.4;y 为 高 度 ;z() 为 y 处 的 风速 值 。 

采用 欧 拉 双 流体 模型 ,附加 应 流 模型 。 气 固 相 
之 间 采 用 Schiller-naumann 计算 相互 间 阻 力 。 由 于 
需要 观测 流 场 随时 间 的 变化 情况 , 故 采 用 非 定 常 的 
瞬 态 求解 方法 来 模拟 ,空间 离散 格式 采用 二 阶 迎 风格 
式 , 时 间 步 长 取 0.001 s, 流 场 求解 采用 SIMPIE 算 法 。 
1.2.3 控制 方程 ”由 于 风沙 流 之 间 热 量 交 换 可 忽略 
不 计 , 故 不 涉及 能 量 方程 。 本 文 模拟 气流 为 不 可 压 
缩 流 , 模 拟 包 含 的 控制 方程 主要 有 连续 方程 .动量 
方程 和 满 流 模型 方程 。 标 准 k-e 模 型 的 稳定 性 、 经 济 
性 和 计算 精度 均 较 高 ,使 用 范围 广泛 ,适合 高 雷诺 
数 的 汕 流 使 用 ,满足 本 文 的 计算 要 求 。 该 模型 需要 
求解 满 动能 及 其 耗 散 率 方程 , 沸 动 能 运输 方程 是 通 
过 精确 的 方程 推导 得 到 ,但 耗 散 率 方程 是 通过 推理 
以 及 数学 上 模拟 相似 原形 方程 得 到 。 汕 动能 大 和 耗 


散 能 se 方程 如 下 : 
Tia AE A ITRE: 
dk _ ð LK, \ ok _ 
Ad 下 人 +o pe (2) 
耗 散 能 se 方程 : 
de _ ð M, \de E€ _ Ee 
P dt i P us tjale, (G,+C,,6,) CD k 
(3) 
eee pC, k 、 Een = 
RP RRE /= 人 全 ，C ,为 经 验 常数 ;G, 表 示 


E 
由 于 水 平 速 度 梯度 引起 的 滑动 能 ;表示 由 于 浮力 
影响 引起 的 满 动 能 ;1 为 时 间 ;xiww% 分 别 表 示 x、y 方 癌 
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的 微 元 体 分 量 ; C,, = 1.44 ; C,,=1.92 ; C,,=0.09 ; 
o,=1.0;0,=13. 

124 防风 效率 ”通过 防风 效率 来 反映 植株 对 气流 
的 削减 作用 ,公式 如 下 : 


na=(1- 2} 100% (4) 


式 中 :x 为 距 植 株 的 水 平 距离 ;z 为 距 地 面 高 度 ; 7 。 
为 (x,z) 处 的 防风 效率 ; v ,为 气流 经 植株 后 在 (x,z) 
处 的 风速 值 ; v, 为 未 经 植株 在 z 高 度 处 的 初始 风 
速 值 。 
1.3 野外 现场 试验 设计 

为 了 验证 模拟 的 可 靠 性 ,在 中 国 科 学 院 策 勒 沙 
漠 研 究 站 实验 观测 场 开 展 野 外 试验 。 试 验 时 间 为 
3 一 5 月 ,此 时 该 地 区 气温 回升 ,风沙 活动 尤为 频 
繁 。 选 择 零 散 分 布 自然 生长 干枯 骆驼 刺 的 沙 地 以 
及 无 植被 生长 的 裸 沙 地 布设 样 地 , 样 地 大 小 为 6 mx 
6m, 其 中 骆驼 刺 样 地 以 骆驼 刺 为 中 心 建 立 ,并 保证 
该 样 地 内 有 且 只 有 2 个 治 主 风 回 分 布 的 间距 为 1 m 
的 植株 。 在 两 植株 中 间 位 置 以 及 裸 沙 地 样 地 中 心 
布设 集 沙 仪 (图 3)。 集 沙 仪 设置 16 个 集 沙 盒 ,每 个 
集 沙 盒 口 径 为 1.5 cmx3 ecm, 收集 距 地 表 0~48 cm 高 
度 内 的 沙 粒 。 集 沙 仪 的 开口 与 观测 期 间 风 向 正 对 ， 
其 底部 与 地 面 齐 平 。 集 沙 时 长 为 24h, 开 始 集 沙 时 
打开 集 沙 仪 进 沙 口 , 集 沙 结束 后 将 其 关闭 并 将 收集 
到 的 沙 粒 融 回 实验 室 , 用 精度 天 平 进行 分 层 称 重 。 

输 沙 率 表示 在 一 定 风 速 和 沙 源 条 件 下 ,地 表 的 
输 沙 能 力 ”。 为 表述 统一 ,本 研究 使 用 输 沙 率 作为 
描述 单位 时 间 .单位 面积 风蚀 物质 量 的 物理 量 , 单 
位 采用 gm”min 。 假 设 某 个 高 度 处 的 集 沙 盒 在 某 
一 时 间 段 7(min) 内 收集 的 沙子 质量 为 m(g) , 则 输 
WKH: 


On. (5) 
式 中 ;5 为 集 沙 合 入 口 面积, 取 值 为 0.00045 m?, 


2 结果 与 分 析 


2.1 流 场 模 型 合理 性 验证 

俞 明 聪 涵 的 研究 表明 ,风沙 运动 中 沙 粒 体积 分 
数 与 输 沙 率 的 分 布 情 况 一 致 。 本 文 利用 野外 试验 
得 到 的 输 沙 率 数据 作为 对 照 ,通过 数值 模拟 获取 植 
PR 1 E AVA 0.5 mm 处 沙 粒 体积 分 数 随 高 度 的 变化 规 
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图 3 集 沙 仪 
Fig. 3 Sand sampler 


律 ,两 者 进行 对 比 验 证 。 图 4、 图 5 显示 ,植株 后 输 沙 
率 和 沙 粒 体积 分 数 随 高 度 的 分 布 基本 一 致 ,在 15~ 
24 cm 高度 处 两 者 均 出 现 波 动 ,呈现 出 先 增 大 后 减 
小 的 规律 ,形成 此 现象 的 原因 与 植株 在 不 同 高 度 处 
的 防风 效率 不 同 有 关 , 同 时 ,植株 后 的 涡流 效应 可 


45 


Ne 一 裸 沙 地 
目 一 - 植株 覆盖 地 表 
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x 15 
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© 5 

0 
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高 度 /m 
图 4 输 沙 率 随 高 度 变 化 (野外 试验 ) 
Fig. 4 Variation of sediment transport rate 
with height (field test) 
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图 5 沙 粒 体积 分 数 随 高 度 变 化 (数值 模拟 ) 
Fig. 5 Variation of sand volume fraction with height 


(numerical simulation) 


能 会 产生 一 定 的 影响 。 就 野外 试验 得 到 的 输 沙 率 
而 言 , 相 较 于 裸 沙 地 近 地 表 处 的 输 沙 率 OE BE 
0.5 mAh 30 cm 高 度 以 下 的 输 沙 率 大 大 降低 ,其 阻 沙 
能 力 良 好 。 验 证 结果 总 体 呈 现 较 好 的 一 致 性 ,因此 
数值 模拟 中 流 场 设置 相对 合理 可 靠 , 后 续 将 采用 相 
同 的 流 场 设置 进行 模拟 。 
2.2 干枯 骆驼 刺 植 株 周围 流 场 特征 

为 研究 植株 影响 下 的 流 场 特征 ,选择 高 度 2 m 
处 风速 为 6m*s!' 的 风速 廊 线 进行 模拟 (下 文 速度 未 
作 特 别 说 明 均 为 2 m 高 度 处 的 风速 ), 流 场 变化 云图 
如 图 6 所 示 。 从 图 6 可 以 看 到 ,植株 明显 改变 了 流 
场 的 分 布 ,气流 形成 分 区 ,为 遇 阻 减速 区 .抬升 加 速 
区 、 亲 流 减 速 区 ,恢复 区 。 在 植株 1 前 1 m 范 围 内 ， 
受 植株 的 阻碍 作用 ,气流 随 靠近 植株 速度 逐渐 减 
小 ;植株 1 处 ,气流 产生 汇聚 并 受到 抬升 ,植株 1 上 
方 风 速 增 大 ,50 cm 高度 处 的 风速 由 4.63 mes 
至 5.01 ms ', 增 大 8.21%, 而 植株 2 在 植株 1 的 有 效 
防护 距离 内 ,其 上 部 并 没有 形成 明显 的 加 速 区 ; 植 
株 1.2 之 间 以 及 植株 2 之 后 一 段 距离 内 风速 减 小 ， 
形成 明显 大 范围 的 亲 流 减速 区 ,该 区 域内 沙 粒 所 获 
动能 减 小 , 易 形成 沉积 ; 随 着 气流 逐渐 远离 植株 区 
域 ,气流 受到 的 扰动 作用 减弱 ,速度 逐渐 恢复 。 

气流 在 植株 处 受到 挤 压 ,越过 植株 后 ,产生 分 
离 ,过 流 断 面 增 大 ,由 于 能 量 差异 ,形成 向 下 运动 的 
分 流 ,靠近 地 面 处 的 气流 返回 入 口 方向 ,从 而 形成 
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遇 阴 减速 区 


0.00e+00 1.04e+00 2.07e+00 


3.11e+00 
风速 /m.sD 


4.14e+00 5.18e+00 6.21e+00 6.90e+00 


图 6 植株 周围 风速 云图 
Fig. 6 Cloud chart of wind speed around the plant 


了 回流 , 即 在 植株 背风 侧 减 速 区 产生 反问 气流 , 形 
成 涡流 现象 ,但 由 于 植株 具有 良好 的 药 透 度 ,该 现 
象 并 不 是 特别 明显 。 由 图 7 可 知 , 流 场 内 形成 3 个 
主要 涡流 区 域 :涡流 1 处 于 两 植株 之 间 的 区 域 , 涡流 
2 处 于 植株 2 至 距 和 人口 7.8 m 处 ,涡流 3 处 于 7.8 m 至 


7.70e-07 1.50e—01 3.00e-01 


涡流 2 


距 入 口 10.5 m 人 处 ,3 个 回流 范围 均 较 大 ,而 强度 均 较 
小 。 涡 流 3 产 生 的 原因 可 能 是 受 涡流 2 末端 的 影 
响 ,气流 再 次 受到 挤 压 后 形成 了 微弱 的 涡流 ,其 后 
气流 逐步 恢复 。 在 涡流 区 域 , 沙 粒 会 被 分 成 两 部 
分 ,一 部 分 沙 粒 随 气流 回流 并 逐步 沉积 , 另 一 部 分 


4.50e-01 
风速 /mm.s-) 


0.00e-01 7.50e-01 9.00e-01 = 1.00e+00 


图 7 植株 周围 风速 矢量 图 


Fig.7 Vector diagram of wind speed around the plant 
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随 涡 旋 气流 进入 气流 主流 区 离开 。 回 流 区 的 范围 
越 大 强度 越 弱 , 越 多 的 沙 粒 会 沉积 在 植株 附近 , 阻 
沙 能 力 越 强 ,这 也 反映 出 植株 在 防风 阻 沙 方面 远 优 
于 无 孔 际 度 的 障碍 物 。 

2.3 干枯 骆驼 刺 植 株 的 防风 效果 

2.3.1 速度 分 布 特征 以 6m's -为 人口 风速 ,风速 
方向 为 正 向 ,提取 0.1m .0.2 m,0.3 m,0.4 m,0.5 m, 
1.5 m 高 度 处 风速 的 水 平 速度 分 量 绘制 沿 程 分 布 图 ， 
其 中 ,植株 实体 位 置 处 风速 取 值 为 0。 由 图 8 可 知 ， 
1.5 m 高 度 处 的 水 平 风 速 一 直 处 于 稳定 状态 ,其 他 5 
个 高 度 处 的 气流 受到 了 不 同 程度 的 扰动 。0.4 m、 
0.5 m 处 的 水 平 风速 在 植株 前 先 略 微 减 小 ,经 植株 失 
升 速 度 增 大 后 ,越过 植株 后 平稳 减 小 ,最 终 在 距 人 
口 14 m 后 逐渐 恢复 到 初始 风速 。0.4 m 人 处 水 平 风 速 
沿 程 最 大 增加 7.87% ,最 大 可 以 降低 39.93%,0.5 m 
处 则 分 别 为 8.29% 17.99% ,增幅 接近 ,降幅 则 呈现 
越 靠近 地 面 越 大 的 规律 。 在 等 于 或 低 于 植株 的 高 


N vU A U A N 


IKE WE/(m-s) 


度 处 (0.1 m、0.2 m、0.3 m) ,气流 受到 的 扰动 增加 ,水 
平 风 速 在 越过 植株 后 呈现 出 不 同 的 变化 规律 。 首 
先 , 不 存在 明显 高 于 初始 水 平 风速 的 阶段 , 即 风 速 
沿 程 整体 呈现 出 被 削减 的 状态 ,这 与 在 低 于 植株 的 
高 度 内 产生 了 有 效 的 防护 有 直接 关系 。 当 风速 减 
小 至 负 值 时 ,表明 此 处 产生 回流 ,附近 范围 内 会 有 
涡流 形成 ,0.1 mm 高 度 处 的 气流 在 较 大 的 范围 内 存在 
回流 现象 ,该 高 度 附 近 是 涡流 现象 的 主要 回流 区 , 
沙 粒 的 沉降 也 主要 在 该 高 度 内 发 生 。 其 次 ,气流 2 
次 越过 植株 后 ,水平 风速 均 呈 现 “W” 型 分 布 ,在 两 植 
株 之 间 ,风速 运行 较为 稳定 。 总 体 而 言 , 除 植株 处 ， 
各 高 度 水 平 风 速 沿 程 变化 整体 较 平 组 ,无 剧烈 扰 
动 , 同 时 在 较 大 范围 内 低 于 入 口 风速 ,有 利于 抑制 
风沙 流 ,促进 沙 粒 沉降 。 

以 6m':s 为 人 口 风 速 ,在 距 入 口 4.3m、S.6m、6 
m,7.2 m,7.8 m、10m 处 取 垂 直 截 面 ,绘制 水 平 风速 
分 量 随 高 度 分 布 图 。 从 图 9 可 以 看 到 ,4.5 mm 处 的 水 


一 o m 


T 


3.00 4.03 4.67 5.00 5.23 5.47 5.79 6.11 6.36 6.59 6.89 7.36 8.00 8.88 10.10 11.76 14.00 
距 入 口水 平 距 离 /m 


图 8 不 同 高 度 气流 水 平 速度 分 量 沿 程 分 布 


Fig. 8 Distribution of horizontal airflow velocity components at different heights 
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图 9 不同 高 度 气流 水 平 速度 分 量 沿 高 分 布 


Fig.9 Distribution of horizontal airflow velocity components along different heights 
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平 风速 基本 以 较为 平稳 的 趋势 增加 ,该 位 置 处 于 植 
株 1 前 0.5 m 处 ,气流 受 植株 影响 相对 较 小 ,变化 的 
幅度 也 相应 的 较 小 。5.6 m、6 m 处 于 两 植株 之 间 ， 
7.2 m 处 于 植株 2 后 0.6m 处 ,3 个 位 置 处 的 水 平 风速 
在 0.6 m 高 度 以 下 均 呈 “W” 型 分 布 , 随 着 高 度 的 不 断 
增加 ,水 平 风速 不 断 减 小 至 负 值 , 即 出 现 回流 , 负 值 
主要 存在 于 0.08~0.14 m 高 度 内 ,并 在 0.12m 附 近 出 
现 第 1 个 极 小 值 , 即 5.6 m、6 m、7.2 m 人 处 分 别 为 : 
-3.3 ms -0.23 m's” -0.15 ms '; 随 着 高 度 增加 ， 
水 平 风 速 增加 后 再 次 减 小 ,并 在 0.27 m 附 近 达 到 第 
2 个 极 小 值 (1.09 ms 1'、1.39 m+s',0.58 mes!) ,该 处 
并 未 产生 回流 。10 m 处 的 水 平 风 速 整体 呈现 随 高 
度 增加 而 增加 的 规律 ,以 0.1 m 高 度 为 节点 ,其 下 为 
回流 ,其 上 为 正 向 流动 ,这 正 是 上 文 分 析 的 较 大 范 
围 的 涡流 3。 总 体 而 言 ,0.6 m 高 度 是 风速 发 展 的 重 
要 分 界 点 ,其 下 风速 先 被 减 小 ,到 达 0.3~0.6 m 人 处 再 
以 较 大 的 加 速度 发 展 ,其 上 风速 几乎 不 受 植株 的 影 
啊 ,风速 恢复 对 数 分 布 , 各 位 置 风 速达 到 一 致 。 
此 ,植株 的 防护 高 度 为 0.6 m, 同 时 0.3~0.6 m 内 风速 
急剧 加 速 的 现象 值得 引起 注意 。 

2.3.2 植株 的 防风 效率 ”风沙 流 的 运动 过 程 中 沙 粒 
主要 集中 在 近 地 面 0.3 m 高 度 内 ,因此 选取 风速 为 6 
ms 10 m's” 时 ,0.1m、0.2m、0.3 m 3 个 高 度 处 的 
防风 效率 进行 研究 。 图 10 显示 ,不 同 风 速 下 防风 效 
率 整 体 呈 现 出 随 着 远离 植株 下 降 的 趋势 。 当 人 口 
风速 为 6m:s"' 时 , 株 后 5.3 m 范 围 内 ,0.3 m 高度 以 下 
的 防风 效率 均 高 于 40% ,各 高 度 处 风速 在 距离 植株 
7.4 m 处 接近 ,为 30% 左 右 ;0.1 m 高 度 处 的 防风 效率 


防风 效率 /% 
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整体 较 高 ,尤其 是 在 距 植 株 3.9 m 范 围 内 ,总 体高 于 
80% ,在 植株 后 1.2~1.6 m、2.8~3.2 m 人 处 出 现 波动 ， 
5.5 m 后 降低 到 50% 以 下 ;0.2 mm 高度 处 的 防风 效率 
在 1.3 m 范 围 内 先 减 小 后 增 大 ,1.3 m 后 防风 效率 逐 
渐 减 小 并 在 1.75 m 人 处 开始 低 于 80% ,植株 后 5 m 降 
低 到 50%。 相 较 于 0.1 m、0.2 m 人 处 ,0.3 m 高 度 处 的 
防风 效率 呈现 出 随 远 离 植株 更 为 稳定 的 减 小 趋势 ， 
防风 效率 低 于 80% ,并 在 2.7mm 处 降低 到 50% 以 下 。 
当 入 口 风 速 为 10 mss 时 ,0.1m 高 度 内 的 防风 效率 
在 近 上 距离 内 良好 ,但 4.26mm 后 迅速 减 小 ,0.2 m .0.3 m 
处 也 有 类 似 的 变化 规律 ; 株 后 5.3 m 范 围 内 ,0.2 m 高 
度 以 下 的 防风 效率 均 高 于 30% ,0.3 mm 处 减 小 到 
16.56%。 总 体 而 言 , 风速 增 大 ,植株 的 防风 效率 减 
小 , 且 该 现象 随 高 度 增加 愈加 明显 ; 当 风 速 较 小 时 ， 
植株 的 防护 范围 较 大 ,防风 效率 较 高 , 当 风 速 较 大 
时 ,贴近 地 面 处 防风 效率 较 高 ,但 在 达到 某 一 距离 
后 其 减 小 速率 增 大 ,植株 较 快 的 失去 有 效 的 防护 ， 
防护 范围 也 相应 减 小 。 
2.4 干枯 骆驼 刺 植 株 周 围 积 沙特 征 

风 是 沙 粒 运动 的 直接 动力 ,风速 的 改变 直接 影 
响 沙 粒 的 运动 状态 。 受 植株 的 影响 ,一 方面 风速 减 
小 , 沙 粒 沉降 , 男 一 方面 ,运动 的 沙 粒 遇 植株 受阻 ， 
形成 积 沙 。 为 研究 不 同 风速 条 件 下 植株 对 积 沙 的 
影响 , 取 入 口 风速 为 6m':s 10m':s 20m:s ,模拟 
44 Bl) =4 s、t=10 s 不 同 风速 条 件 下 的 积 沙 云图 (图 
11)。 从 图 11 可 以 看 到 ,不 同 风速 条 件 下 , 随 着 运动 
时 间 的 增加 , 积 沙 会 不 断 增 加 ,并 在 某 一 时 刻 达 到 
该 风速 下 的 稳定 状态 。 当 入 口 风速 为 6m's Ht, Yb 


—— 0.1 m (6 ms’) 
一 一 0.2 m (6 ms") 
一 一 0.3 m (6 m-s”) 
—— 0.1m (10 ms") 
goes —— 0.2 m (10 ms") 
一 一 0.3 m (10 ms") 


0.10 0.29 0.50 0.76 1.05 140 1.81 2.28 2.84 3.50 426 5.16 6.22 7.40 


与 植株 2 水 平 距离 /m 


10 防风 效益 
Fig. 10 Protection benefit 
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(a) f=4 s, v=6 mes! 


(c) f=4 s, v=10 mes 


(e) f=4 s, v=20 mes! 


0.00e+00 1.04e+00 2.07e+00 


3.11e+00 
风速 /(m.s-!) 


(b) 10 s, v=6 mes? 


4.14e+00 5.18e+00 6.21e+00 6.90e+00 


注 : 不 同 的 体积 分 数 通过 颜色 区 分 ;红色 为 最 大 值 , 表 示 沙 粒 在 该 处 有 沉积 现象 ; 蓝 色 为 最 小 值 ; 
其 他 颜色 代表 沙 粒 以 一 定 的 形式 (里 移 RAS RAS HEFT IZ 
图 11 RA BARU Zs Al 


Fig. 11 Sand accumulation cloud atlas around the plant 


粒 先 在 植株 1 附近 、 植 株 2 迎 风 侧 处 沉积 ,随后 在 植 
株 1 处 及 两 植株 之 间 不 断 积累 ,植株 1 附近 的 积 沙 
明显 大 于 植株 2 附近 的 积 沙 ,这 是 由 于 气流 遇 植株 1 
受阻 .风速 降低 .搬运 沙 粒 能 力 下 降 导 致 的 ,与 前 文 
风速 分 析 结 果 相 对 应 。 当 入 口 风速 为 10 m's 时 ， 
同样 是 在 植株 1 附近 及 两 植株 之 间 形 成 积 沙 ,但 当 
t=10 s 时 ,植株 1 附近 的 积 沙 相 较 于 同时 刻 入 口 风速 
H 6 m's 时 减少 ,更 多 的 积 沙 出 现在 植株 2 前 。 随 
着 入 口 风速 增加 到 20 m:'s ,在 2 个 时 间 点 均 呈 现 植 
株 2 处 形成 更 多 的 积 沙 的 现象 ,尤其 是 当 大 10 s 时 ， 
植株 1 附近 的 积 沙 与 /=4 s 时 无 明显 差别 , 均 存 在 少 
量 积 沙 , 而 植株 2 前 后 形成 大 量 积 沙 , 其 至 形成 小 沙 
堆 。 分 析 其 原因 如 下 :风沙 流 经 过 植株 时 气流 速度 
发 生 明 显 改 变 , 当 风速 较 小 时 ,植株 1 的 削减 作用 ， 
使 其 前 后 一 定 范围 内 的 流速 降低 到 起 沙 风 速 以 下 ， 
风沙 流 处 于 过 饱和 状态 , 沙 粒 沉降 堆积 ,但 当 风 速 
较 大 时 ,虽然 植株 1 对 风速 有 削减 ,但 减弱 能 力 有 
限 ,风速 仍 大 于 起 沙 风速 ,风沙 流 处 于 非 他 和 状态 ， 
仅 小 部 分 沙 粒 受 阻 沉积 ,大 部 分 沙 粒 继续 运动 , 植 
株 2 对 气流 形成 二 次 削减 ,附近 风速 大 大 降低 ,此 时 
大 量 沙 粒 沉积 ,形成 积 沙 ,同时 ,植株 2 后 较 大 范围 
的 回流 也 促进 了 沙 粒 的 沉积 。 

研究 显示 ,植株 具有 良好 的 阻 沙 能 力 。 在 野 


外 , 随 着 时 间 的 推移 ,不 同 速度 风沙 流 交替 运行 , 植 
株 附近 沙 粒 不 断 堆 积 ,植物 也 通过 生理 作用 适应 风 
沙 活 动 ,经 过 长 时 间 的 相互 作用 ,植株 附近 形成 沙 
堆 , 也 就 是 常 说 的 灌 从 沙 堆 。 灌 从 沙 堆 或 大 或 小 ， 
从 几 十 厘米 到 几米 不 等 ,这 与 植株 的 大 小 、 风 沙 活 
动 的 强 弱 、 当 地 环境 特征 等 多 种 因 系 有 关 。 由 于 野 
外 风速 多 变 , 较 小 风速 下 ,前 面 的 植株 可 使 沙 粒 沉 
降 , 但 当 风 速 较 大 时 , 则 需要 后 面 的 植株 进行 全 加 
作用 ,因此 多 植株 比 单 植株 更 易 形成 灌 从 沙 堆 , 且 
其 会 受到 植株 布局 的 影响 。 通 过 栽种 植物 或 者 利 
用 类 似 于 植株 的 材料 进行 防风 治 沙 时 ,要 考虑 到 区 
域 的 风速 特征 ,在 易 形 成 高 风速 风沙 流 的 地 区 ,应 
尽量 布设 双 行 甚至 多 行 来 达到 较 好 的 效果 。 


3 讨论 


引起 风沙 灾害 .影响 土壤 风蚀 强度 的 因素 有 很 
多 ,包括 风 场 .土壤 性 质 . 地 表 条 件 等 ,植被 覆盖 是 
地 表 条 件 的 主要 影响 因子 之 一 。 由 于 区 域内 风 场 、 
土壤 性 质 具 有 固定 性 ,利用 植被 改变 地 表 条 件 以 达 
到 防风 固沙 、 抑 制 风蚀 的 目的 是 十 分 有 效 的 措施 ， 
因此 对 植被 影响 下 流 场 特征 的 研究 尤为 关键 。 但 
是 ,由 于 植被 及 其 影响 下 的 地 表 极 其 复杂 ,不 同 覆 
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TEE Ty SE PR ee eat a PE .配置 模式 的 植被 对 地 表 物 
质 条 件 及 风速 的 影响 均 有 差异 ”, 使 得 该 研究 十 分 
困难 。 

植被 主要 通过 2 个 机 制 对 风沙 活动 产生 影响 ， 
学 者 们 根据 侧重 点 不 同 ,针对 这 2 个 方面 分 别 展开 
人 研究。 一 方面 ,植物 改变 风 场 特征 ,降低 风速 ,同时 
对 沙 粒 造 成 阻碍 , 沙 粒 沉降 。 徐 高 兴 等 ”通过 风 洞 
试验 ,模拟 了 4 种 配置 模式 下 梭 梭 林 (Haloxylon am- 
modendron ) 的 固沙 效果 ,人 研究 表明 不 同 配 置 模式 会 
影响 流 场 运 行 特征 ,改变 固沙 效果 ,但 当 植 被 覆盖 
度 达到 32.37% 时 ,各 配置 模式 下 梭 梭 林 的 固沙 效果 
差异 不 显著 。 对 于 小 型 灌木 , 受 植株 形态 的 影响 ， 
术 下 高 度 大 的 植株 形成 的 " 狭 管 效 应 ”会 加 速 地 表 的 
风蚀 作用 ,形态 低 矮 且 玻 透 度 小 的 植株 防风 阻 沙 效 果 
更 明显 ,有 效 防 护 距 离 更 大 ” 。 另 一 方面 ,植物 通过 
生物 作用 使 土壤 形成 有 机 质 层 ,固定 土壤 , 减 小 地 
表 剪 切 力 ,提高 地 表 起 沙 风 速 ,降低 沙 粒 流动 幅度 ， 
从 而 达到 防风 固沙 的 效果 。 王 佳 庭 等 ”对 乌 兰 布 
和 沙漠 不 同 植物 群落 影响 下 ,土壤 可 蚀 性 因子 的 研 
究 发 现 ,不 同 沙 地 可 蚀 性 颗粒 含量 ` 有 机 质 含量 E 
REAT .土壤 含水 率 受 植被 影响 存在 差异 , 梭 
权 AH (Nitraria tangutorum ) 等 灌木 植物 群落 比 沙 
3 (Agriophyllum sguarrosumz)、 盐 爪 爪 人 Kalidiam foli- 
atum ) 等 草本 植物 群落 更 能 降低 土壤 风蚀 作用 ; 吴 
楠 等 的 研究 表明 ,植物 作用 下 土壤 有 机 质 、 养 分 
促进 土壤 结 皮 的 形成 和 发 育 , 使 土壤 具备 更 强 的 蓄 
水 能 力 和 抗 风蚀 能 力 。 数 值 模拟 的 方法 可 以 较 好 
地 反映 出 植株 第 1 个 防风 阻 沙 机 制 的 防护 效果 ,但 
难以 综合 反映 植物 的 生物 作用 。 

不 同 物候 期 ,植株 的 冠 层 形态 ARARE 
面积 以 及 由 此 导致 的 地 上 生物 分 配 模式 具有 差异 ， 
其 阻 沙 能 力 不 尽 相同 ”, 而 灌木 植株 附近 沙 堆 的 高 
度 受 植物 高 度 的 控制 , 沙 堆 的 半径 受 冠 幅 的 控制 ， 
沙 堆 的 生长 发 育 受 形态 参数 的 影响 ”。 基 于 生长 
季 骆 驼 刺 与 非 生 长 季 骆 驼 刺 在 高 度 、 冠 幅 及 孔 院 度 
方面 的 差异 ” ,本 文 的 植株 模型 通过 较 小 的 高 度 和 
冠 幅 、 较 大 的 孔 际 来 表现 非 生长 季 骆 驼 刺 的 形态 有 具 
有 一 定 的 可 靠 性 。 后 续 人 研究 中 需 优化 生长 季 骆 驶 
刺 植株 模型 ,同时 考虑 植物 生物 作用 以 及 野外 非 定 
常 来 流 因素 的 影响 ,以 获得 更 全 面 、 准 确 的 结果 。 
同时 ,模拟 结果 显示 出 2 个 植株 附近 这 一 极 小 尺度 

下 流 场 的 运行 情况 ,大 尺度 背景 下 不 同窗 盖 度 植株 
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的 协同 作用 是 男 一 项 重要 的 人 研究 内 容 “,2 个 尺度 
关系 密切 ,小 尺度 环境 影响 大 尺度 流 场 的 运行 ,大 
尺度 又 可 进一步 反映 小 尺度 ,因此 进行 相关 人 研究 对 
建立 两 者 相互 佐证 十 分 关键 。 


4 结论 


本 文 对 30 cm 高 度 骆 驼 刺 的 形态 进行 简化 , 建 
立 二 维 模型 ,运用 数值 模拟 的 方法 对 风沙 流 经 过 植 
株 后 的 风速 特征 和 积 沙 特征 展开 人 研究 ,并 利用 对 外 
试验 验证 其 可 靠 性 ,研究 结论 如 下 : 

(1) 植株 附近 流 场 大 致 可 以 分 为 遇 阻 减速 区 、 
抬升 加 速 区 . 北 流 减速 区 和 恢复 区 。 植 株 周围 会 形 
成 3 个 微弱 的 涡流 ,分 别 位 于 两 植株 之 间 、 植 株 后 
1.2m 内 以 及 1.2~2.7 m 内 。 人 靠近 植株 处 涡流 的 回流 
区 主要 存在 于 0.09~0.14 m 高度 内 ,远离 植株 处 的 回 
流 区 位 于 0.1m 以 下 。 

(2) 当 风 速 为 6m's AY, 30 cm 高 度 植株 主要 影 
响 0.6mm 以 下 的 水 平 风速 ,在 离开 两 植株 8 mm 左右 恢 
复 到 初始 状态 ,植株 附近 的 风速 沿 程 呈 “W” 型 分 
布 。 植 株 后 一 定 距 离 内 的 水 平 风 速 随 高 度 增加 存 
在 2 次 减 小 过 程 ,并 在 0.3~0.6 m 高 度 范围 内 , 沿 高 
度 以 较 大 的 加 速度 增 大 后 ,逐步 恢复 至 初始 风速 廓 
线 状态 。 

(3) 植株 的 防风 效率 整体 呈现 出 随 风 速 增 大 而 
减 小 的 规律 , 量 该 现象 随 高 度 增加 愈加 明显 。 由 于 
植株 具有 良好 的 玻 透 度 ,其 防风 效率 整体 良好 , 防 
护 距 离 较 大 , 当 风 速 为 6m:s' 时 ,植株 后 5.3 m 范 围 
内 ,0.3 m 以 下 的 防风 效率 高 于 40%, 当 风速 增加 至 
10 ms" 时 ,该 值 减 小 至 16.56%。 

(4) 植株 的 阻 沙 能 力 良 好 。 不 同 风速 条 件 下 ， 
随 着 时 间 的 增加 , 积 沙 会 不 断 增 加 并 达到 稳定 。 当 
风速 较 小 时 , 积 沙 主要 集中 在 前 方 植株 附近 以 及 植 
株 之 间 的 区 域 , 随 着 风速 的 增 大 ,气流 经 多 次 削减 
后 沙 粒 沉积 , 积 沙 有 后 移 的 趋势 。 

致谢 :本 次 野外 实验 的 数据 获取 得 到 了 中 国 科 
学 院 策 勒 沙漠 研究 站 的 大 力 支 持 , 对 此 表示 衷心 的 
感谢 ! 
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Numerical simulation study on the influence of dry Alhagi 
camelorum on the wind-sand flow field 


LIU Jinmiao’, LI Juyan’, YIN Zhongdong', GUAN Hanxiao’, ZHANG Jiawei’ 
(1. School of Soil and Water Conservation, Beijing Forestry University, Beijing 100083, China; 2. Xinjiang General 


Ecological Environment Monitoring Station of Soil and Water Conservation, Urumqi 830002, Xinjiang, China; 
3. School of Life and Geography, Kashi University, Kashi 844000, Xinjiang, China) 


Abstract: One of the important measures to control sand disasters in western desert areas is sand fixation by 
plants. As a typical desert plant, Alhagi camelorum has important application value. Based on the flow field in the 
spring, 30 cm of dry Alhagi camelorum was the object of research in this paper. To analyze the characteristics of 
wind speed and sand deposition, a numerical simulation was used via the fluent and field tests. The results were 
as follows: (1) the wind-sand flow can be roughly divided into the areas of blocked deceleration, lifting acceleration, 
turbulent deceleration, and recovery when going through the plant. A weak vortex forms behind the plant, and the 
height of the recirculation region is related to the distance from the plant. However, they are less than 0.14 m. (2) 
When the wind speed is 6 m*s , plants at a height of 30 cm mainly affect the horizontal wind speed below a 
height of 0.6 m. The horizontal wind speed at a certain distance behind the plant no longer presents a strict 
logarithmic distribution with height, but there are two minima, and the acceleration increases rapidly in the height 
range of 0.3-0.6 m. (3) The windproof efficiency of plants decreases with the increase of wind speed, and this 
phenomenon becomes more and more obvious the height increases. When the wind speed increases from 6 to 10 
m's ,the windproof efficiency at a height of 0.3 m within a distance of 5.3 m behind the plant decreases from 
40% to 16.56%. (4) The sand accumulation near the plant is different due to the different wind speed. When the 
wind speed of the incoming flow is low, sand accumulation is mainly concentrated near the front of the plant and 
within the range between plants. As wind speed increases, sand accumulation moves backward. 


Keywords: sand-fixing service; flow field distribution; numerical simulation; Alhagi camelorum; vortex 


